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Q^ I A number of recent determinations of m^ using hadronic r decay data involve 

[~^ ' inclusive analyses based on the so-called (k, 0) spectral weights. We show here 

^^ ■ that the OPE representations of the longitudinal contributions appearing in 

f— ^ . these analyses, which are already known to be poorly converging, have in 

addition an unphysical k dependence which produces a significant unphysical 

0^1 decrease in nig with increasing k. We also show how, using additional sum 

o ' rule constraints, the decay constants of the excited resonances in the strange 

^ ■ scalar and pseudoscalar channels may be determined, allowing one to evaluate 

• • . the longitudinal spectral contributions to the {k, 0) sum rules. Taking into 

. ^H ! account the very-accurately known vr and K pole contributions, we find that 

kS ' longitudinal contributions can be determined with an accuracy at the few 

C^ ' % level, and hence reliably subtracted, leaving an analysis for m^ involving 

the sum of longitudinal and transverse contributions, for which the OPE 

representation is much better converged. 
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I. INTRODUCTION 

The determination of the strange quark mass, rris, has been the focus of much recent 
activity in both the sum rule [p^-[T^ and lattice |[T5|-|T8| communities. (A summary of the 
current status of both the sum rule and lattice determinations is given in Ref. |]T9| .) 

Among the sum rule approaches, those based on flavor breaking in hadronic r decay [p|-p^ 



appear the most reliable at present. There are two reasons for this statement. First, on the 
experimental side, the spectral data required is known over the full kinematic range entering 
the relevant sum rules ||20|j2T[| and, second, on the theoretical (OPE) side, flavor-independent 
instanton and renormalon effects, which create potential uncertainties in analyses of the 
strange scalar and pseudoscalar channels |l|,@,H, cancel in forming the flavor-breaking r decay 
difference 0. The r decay sum rules are, however, not without complications. The primary 
complication has to do with the behavior of the OPE representation of the contributions 
to the inclusive decay rate of hadronic states with total spin J = 0. In order to elaborate 
on this point, and to fix terminology and notation, we briefly review the relation of rus to 
hadronic r decay. 

The ratio of the inclusive hadronic r decay width through the flavor ij = ud or us vector 
(V) or axial vector (A) weak hadronic current to the corresponding electron decay width. 
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where (7) indicates additional photons or lepton pairs, can be expressed in terms of weighted 
integrals over the spin J = (longitudinal) and J = 1 (transverse) components of the 
corresponding V or A spectral functions |^, where the spectral functions are defined, as 
usual, by 
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In Eq. d^), Ilj^^.y^ are the spin J scalar components of the usual flavor ij V, A current-current 
correlators, 
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Working with the combinations n(°+^)(g^) = U^^^q'^) + Il^^^{q'^) and g^n(°)(g^) which have 
no kinematic singularities, Rij-v,A can then be written p3| , p4| 
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where s = q = —Q, Sew = 1-0194 represents the leading electroweak corrections 25 



Vij are the usual CKM matrix elements, and the second line follows from the first as a 



consequence of Cauchy's theorem. The second hne of Eqs. (H) allows Rij;v,A to be evaluated 
using techniques based on the OPE and perturbative QCD p2|- ^ , ^ ,p7[] . The weights, 



wxiy) = (1 — y) (1 + 2?/) and WL{y) = —2y{l — y) (with y = s/m^) multiplying 11*^°"'"^) 
and Yi^^\ respectively, have (double) zeros at s = m?^^ reflecting the fact that s = m^ lies 
at the edge of hadronic phase space. The resulting suppression of contributions from the 
region of the circular contour near the timelike real axis (the region of potential breakdown 
of the OPE) is responsible for the high quality of the OPE representation of the inclusive 
hadronic rates (see Ref. |^ for a recent review). 



Defining R^ = Rij-y + -Rjj;A, it is then evident that 

Rud R 

^ ud\ \^ us 
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vanishes in the SU{?,)f limit. Defining AH^-^) = Yi^;!} - Yi^^} and Ap(-^) = p^^^ - p'^^) , the 
mass-independent D = contributions, therefore, cancel by construction on the OPE side 
of the sum rule for Ai?^ analogous to Eq. (|) , 
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Neglecting m^ ^ and agrriu^^ms relative to m^, the leading D = 2 term is proportional to ml 
PJ2^. The integrand on the LHS of Eq. (^ is known, as a function of s, from the work of the 



ALEPH collaboration [^,^. On the OPE side, the D = 2 contribution is known in terms 
of ttsi once TUs is fixed; the D = A contribution is known in terms of {rrisss)] and the D = 6 
contribution is small as a result of the cancellation between the dominant D = 6 four-quark 
condensate terms which occurs in the V+A sum [^^. Alternate sum rules which also 
allow the spectral side to be evaluated using the measured s-dependent ud and us number 
distributions, without necessitating a J = 0/ J = 1 separation of the experimental data, can 
be constructed by multiplying the integrands appearing on both sides of Eq. (^) by a common 
analytic factor. For the case that this factor is (1 — y)^y"', the resulting sum rule is said 
to involve the {k,n) spectral weight, and the resulting analogue of ARr is denoted AR['''"-\ 
The {k, 0) spectral weight sum rules (with k = 0,1,2) form the basis of a number of recent 
inclusive treatments of the m^ extraction problem [|7|- [lT| , p!3| , p!4| . We will denote the weights 
accompanying 11'^°"'"^) and n*^") in the (/c, 0) sum rules by Wrp (y) = {1 — y) ^ (1 + 2y) and 
^L \y) = ^2?/ {1 — y) ^ , respectively. 

Eq. (|^) (corresponding to the (0,0) spectral weight), and/or its (fc, 0) generalizations, 
would allow a straightforward determination of irts, provided the OPE representations of 
both the J = + 1 and J = contributions above were well converged at scale rn^. Unfor- 
tunately, it turns out that this is not the case. The problem lies with the OPE contribution 
involving the product w}^ {y)I\''^\s). We refer to these contributions as "longitudinal" in 
what follows. 

The source of the problem is that the perturbative series for the integrated D = 2 
longitudinal contribution (which is known to four loops, i.e., 0{a^) p8| , p9| ) is not convergent 
at the scale m^ P,[7|,|30| . This is true whether one considers the "fixed order" (FOPT) 



expansion (expansion in powers of as(At^) at a fixed scale, /i, e.g., /i = m^), or the "contour- 
improved" expansion (CIPT) ^,311 (in wliicli tlie large logarithms are summed up by the 
scale choice /i^ = Q^ point-by-point over the contour). Taking the unmodified version of 
Eq. (H), corresponding to fc = 0, to be specific, one finds that, using the central value of the 
ALEPH determination, as{m1) = .334, the FOPT D = 2 series behaves as 

~ 1 + 0.99 + 1.24 + 1.59 + ■■■ , (7) 

while the corresponding CIPT series behaves as [^ 

~ 1 + 0.78 + 0.78 + 0.90 + • • • (8) 

(where in both cases we have normalized successive terms in the integrated series to the 
leading, 0{a^), term). The integrated longitudinal D = 2 series, truncated at 0(af), also 
exhibits a very strong residual scale dependence 0. Because of the non-convergence and 
strong residual scale dependence, inclusive sum rules of the type described above contain sig- 
nificant uncertainties associated with the presence of the longitudinal contributions. Recent 
inclusive analyses [0,^|-[TT|,T^,n,2T| proceed by taking the sum oiD = 2J = + l and J = 



contributions to O(a^). The size of the last known (0(a^)) term is taken as an estimate of 
the D = 2 truncation error. In Refs. [[7|, [l0| , |l3i p!4| , an additional error has been included to 



account for the residual scale dependence of the truncated sum. The scale-dependence error 
is estimated by varying the scale choice in the CIPT evaluation of the D = 2 sum according 
to /i^ = C,^Q^, with .75 < ^ < 2, and symmetrizing the resulting variation about the central 
value C, = ^- The resulting estimated D = 2 error is much larger than that for the remaining 
longitudinal terms, and hence dominates the error on the total longitudinal contribution. 
Taking the results of Ref. [1^, which gives the detailed breakdown into individual longitu- 



dinal contributions, to be specific, and combining the quoted errors in quadrature, one finds 
errors of 32%, 34% and 37% for the (0, 0), (1, 0) and (2, 0) total longitudinal contributions, 
respectively. These errors represent the dominant component of the total theoretical error 
in the inclusive analyses. 

While the errors on the longitudinal contributions discussed above might appear safely 
conservative, we will argue in this paper that they are, in fact, almost certainly too small. 
We will, in addition, demonstrate that the central values of the OPE contribution to AR^^'^"^ 
contain an unphysical dependence on k which produces a corresponding unphysical lowering 
of the extracted value of rris with increasing k. Finally, we will demonstrate that it is possible 
to significantly reduce the uncertainties on the total longitudinal contributions, allowing one 
to subtract the longitudinal contributions from the experimental number distribution and 
work instead with sum rules for the much better behaved + 1 correlator difference. 

The rest of the paper is organized as follows. In Section II we demonstrate the existence 
of the problem and investigate its magnitude. In Section III we discuss how to improve the 
estimate for the longitudinal contributions using sum rules for the strange scalar and pseu- 
doscalar channels and present our numerical results. Finally, in Section IV, we summarize 
our results and comment on their implications for future analyses of rris using r decay data. 



II. PHYSICAL CONSTRAINTS ON THE LONGITUDINAL CONTRIBUTIONS 
TO THE SPECTRAL WEIGHT ANALYSES 



If a J = 0/J = 1 spin decomposition existed for the current experimental data, one 
could simply subtract the longitudinal contribution from each bin of the experimental ud 
and us distributions, determine Pudulv+Ai^)^ ^^^ use this information to analyze the + 1, 
rather than inclusive, sum rules. Unfortunately, such a decomposition does not yet exist over 
the whole of the kinematically allowed range. Certain general features of the longitudinal 
contributions to the V and A correlators, however, allow us to, nonetheless, obtain useful 
constraints. 

In the chiral limit, the longitudinal spectral functions vanish except for the (massless) vr 



.(0) 



(0) 



and K pole terms, which contribute to p^^-a ^^^ Pus-a-> respectively. Away from the chiral 



.(0) 



.(0) 



and {rrii + m 



limit, pi^y and p^^.^ receive additional contributions proportional, respectively, to (rrii — mj) 

For ij = ud, these additional contributions are numerically tiny and can 
be neglected. For ij = us, the thresholds for the non-pole contributions to p^lv ^^^ PuIa 
are {rriK + "^7r)^ = sfj^ and {tjik + ^rriT.y = s^^"^, respectively. The sum of the longitudinal 



K and vr pole contributions for the {k, 0) spectral weight is 
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with /tt = 92.4 MeV and fx = 113.0 MeV ||3^. The vr pole contribution is nearly constant 
with increasing k, the K pole contribution slowly decreasing with k. The remaining longi- 
tudinal contribution, which we will refer to as the "resonance contribution", is then given 
by 
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where yth = s^^ /m^. The V part of this contribution should be dominated by the i^Q (1430) 
resonance, since only the tail of the next strange scalar resonance, the Kq{1950), lies within 
the kinematically allowed range, s < m^. Similarly, the A contribution in Eq. (p!0| ) should be 
dominated by the i^(1460). The K and vr pole longitudinal contributions are, of course, very 
accurately known, so it is the absence of an experimental determination of the i^Q(1430) 
and ^(1460) decay constants which prevents us from performing a reliable longitudinal 
subtraction. 

Note that, because p^ly and p^Ia are positive, and the analogous ud resonance con- 
tributions are negligible, the longitudinal resonance contribution of Eq. (|T^) is necessarily 
positive and a decreasing function oi k. In fact, from Eq. (|TOD it follows that the longitudinal 
{k, 0) resonance contributions must satisfy the inequality 
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where the equality would obtain only if the entire non-pole us spectral strength lay at thresh- 
old in the scalar channel. That both the pole and resonance contributions are decreasing 



functions of k, of course, also means that the total longitudinal contribution to AR^''^'' must 
be a decreasing function of k. 

The fact that one expects only very small contributions from the tails of the i^Q (1950) and 
ii'(1830) resonances, and that both the masses and widths of the i^o(1430) and i^(1460) 
are similar allows us to sharpen considerably the constraint represented by Eq. ([TI|). In 
the narrow width approximation (NWA), dominance by resonance contributions at M ~ 
1.4 GeV would mean that 
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A more refined version of this estimate is obtained by considering i^Q(1430) and i^(1460) 
Breit-Wigner forms with PDG2000 values for the masses and widths, and integrating directly 
over the resonance profiles. The results of this exercise are that the individual resonance 
contributions to \AR^^^^A for fc = 0, 1,2 are in the ratios 1 : 0.46 : 0.24 for the A:*(1430) 
and I : 0.43 : 0.20 for the i^(1460). We would thus expect the total resonance contributions 
to [Ai?(.'='°)l to satisfy 
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In what follows, we will, for 



Let us now consider the OPE representation of 
convenience, quote the OPE results obtained in Ref. [|T^], since that reference provides a 
complete breakdown of the individual + 1 and longitudinal contributions (in addition to 
a detailed discussion of the evaluation of the various contributions, to which the interested 
reader is referred for details). The D 
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where e^ = rrid/m, 

asim"^). The results of Ref. |iy] are 



.053, and Af^QX, which results from the CIPT integration, depends on 
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= 5.3 ±2.5 ±0.7 


^(2,0) 


= 5.8 ±3.2 ±0.8 



(15) 

where, in each case, the first error represents the combination of the scale-dependence and 
truncation errors, discussed above, and the second represents the effect of the experimental 
uncertainty in ^^(m^). The D = 4 longitudinal contributions are, numerically, completely 
dominated by the term proportional to {rrisss) . Since this term arises from the Ward identity, 

?^ ^IfyA^^) = i^i ± ^jf ^ij;sA(l^) + i^i ± ^j) {{mi) ± {mi)) ' (16) 

where IIjj.s p are the correlators of the flavor ij scalar and pseudoscalar densities, and the 
plus (minus) signs on the RHS correspond to the A (V) case, the Wilson coefficient of the 



{msSs) term receives no radiative corrections. The D = 4 contribution can therefore be 
evaluated rather accurately, using (i) the quark mass ratios obtained by Leutwyler [^], (ii) 
the GMOR relation, ((m„ + md)uu) = -/^m^, [] and (iii) {ss)/{uu) = 0.8 ± 0.2 [0. The 
value turns out to be the same for the (0, 0), (1, 0), and (2, 0) spectral weights [|0 



^Ji(k,0) 
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0.0726 ±0.0194 . 
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2 contributions. The 



The error in Eq. (|T^ is considerably smaller than that on the D 

leading four-quark D = 6 contributions are absent from 11^^., 

D = 6 contributions, while contributions oi D = 8 and higher are assumed to be negligible. 



24], justifying the neglect of 



^Rikfi) 



employed in recent inclusive analyses thus consists 



The OPE representation of 
of the sum oi D = 2 and D 
D = 2 series. 

It is now straightforward to see that the OPE representation just described for 
rather badly violates the constraints obtained above. We first observe that, from Eqs. (p!4D 



4 terms, with the error dominated by that on the truncated 
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and (|T^), AR^'''^^ is a slowly increasing function of k. Since AR^'''^'^ 
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is constant 
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OPE 



is also slowly increasing. 



with k, this means that the OPE representation, 

rather than decreasing, with k!^ Since, for values of ms{rn^) typical of those obtained in 
recent sum rule analyses, ms{m'^) ~ 120 MeV, the D = 2 contributions are more than 
a factor of 2 larger than the corresponding D = 4 contribution, this problem is likely to 
have important numerical consequences. The increase with k of the central values shown 
in Eq. (0), in fact, means that the extracted values of rris must necessarily display an 
unphysical decrease with k. f\ This unphysical decrease would not, of course, be a practical 

(as opposed to conceptual) difficulty for the inclusive analysis if AR). '" represented 

only a small fraction of AR^-. Unfortunately, this is not the case. To illustrate this point 
we show, in Table |, for each of the three {k, 0) spectral weights, the numerical values of the 
longitudinal D = 2 contribution and total longitudinal contribution obtained in Ref. [10 



together with the 1999 ALEPH experimental values of Ai?![.^'°^ [^. The reader should 
be aware that significantly different central values of ms{m'^) were obtained for the three 
different spectral weights; the tabulated D = 2 contributions correspond to these central 
values (also listed in the table). The quoted D = 2 errors were obtained by combining the 



^Deviations from the GMOR relation have recently been shown to be at most 6% |34]. The 
resulting error on the nis analysis is completely negligible. 
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^The small 0{7n'^) contributions to 
also increasing with k, so the full longitudinal D 
increasing with k. 



which have been neglected above, are actually 
= 4 contribution is actually also (very slowly) 



^It is worth noting that Ref. [Q| employs a fc-dependent truncation scheme, in contrast to the 
other inclusive analyses, which truncate at the same order for all k. The problem of the unphysical 
fc-dependence of the central values outlined above, however, remains present even for this altered 
scheme. 



two theoretical errors in Eq. ([ISD in quadrature. We see from the table that the longitudinal 
contribution, in each case, represents more than half of the experimental value. 

Let us attempt to quantify how large the errors associated with the unphysical k depen- 
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tudinal vr and K pole contributions, the OPE representation of 

of as providing an estimate for the sum of the unknown longitudinal resonance contributions. 

r -, ,-| resonance 

In Table || we display the values of AR). ' ' implicit in the OPE results of Ref. []10 

together with the values of the corresponding pole contributions, 
of the resonance contributions are given. In the first (labelled "IIN 
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D" in the table), the values 
oi ms^rn'^) used to evaluate the three different AR^'^'^ are different, corresponding to 

the central values obtained in the relevant independent {k,0) analysis of Ref. |T^. In the 
second (labelled "COMB" in the table), the common value, rris 
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119 MeV, obtained 
T0| , is used to compute all three of the 

with fc, combined with the decrease 
of the sum of the vr and K pole contributions with k, means that the nominal resonance 

OPE 



in the combined fit to all three (fc, 0) sum rules 
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contribution implicit in 
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must be increasing with /c, for fixed ms 



m: 



This is 

evident in the results of the combined analysis, but obscured by the decrease of msijnl) with 
k for the independent analysis. It is evident that both sets of results are far from satisfying 
the constraints given in Eqs. ([T3|) . 

We are now in a position to illustrate the potential significance of the unphysical k 



dependence of 
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OPE 



on the extracted values of mg. Let us imagine that the central 
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and use 



OPE value provides a good approximation for one of the three 
Eq. ( P^ to estimate the resonance contributions to the other two (/c, 0) sum rules. We 
find that, if we attempt to make this assumption for either k = 1 oi k = 2 the result is 
a fc = longitudinal contribution which exceeds the full experimental value, violating the 
positivity of the + 1 OPE representation. If we instead make the assumption for A; = 0, 
the resulting change in the (1, 0) and (2, 0) longitudinal contributions produces a shift in the 
extracted central ms{m'^) values from 121 —>■ 142 MeV and 106 — > 133 MeV, respectively. 
We stress that this exercise is for illustrative purposes only; although the consistency of 
the three analyses is significantly improved if one assumes that, for some reason, the k = 
representation is good, the resulting "extraction" of m^ is meaningless since the assumption 
simultaneously forces the k = 1,2 representations to be bad, leading one to the conclusion 
that it is unreasonable to have assumed that the k = representation was good in the first 
place. 

At present there is little experimental information available on the size of the longi- 
tudinal resonance contributions. The PDG2000 compilation provides no information on 
r —>■ 7^(1460)//,-, and quotes only an upper bound, B < .0005 on the r — >■ i^Q (1430)z/^ branch- 



Ai?(0'0) 



A-*(1430) 



< 0.052. 



ing fraction. The latter bound corresponds to an upper bound 
The central longitudinal (0, 0) OPE determination, if reliable, would then require a corre- 
sponding iir(1460) contribution greater than ~ 0.10. Taking a Breit-Wigner i^(1460) form 
with PDG2000 values for the mass and width, this corresponds to fKiueo) > 100 MeV. 
Such a large value is extremely unnatural given that fK{iim)/ Ik ^ in the SU{?>)p chiral 



r en n^^ -^(1460) 
limit. Not surprisingly, therefore, the lower bound AR)^''' > 0.10 turns out to be 

L -i 1j 

more than an order of magnitude larger than the value obtained from the sum rule analysis 
of the next section. As such, it is completely incompatible with the sum rules for the us 
pseudoscalar correlator. 



III. THE EXCITED STRANGE SCALAR AND PSEUDOSCALAR RESONANCE 
DECAY CONSTANTS FROM SCALAR AND PSEUDOSCALAR SUM RULES 

From Eqs. ( JIB] ) and the Ward identities for the divergences of the flavor ij V and A 
currents, it follows that 

(fpf^-A^i^) = i^i + ^j? PirA<f) = Pir,dA{q^) , (18) 

where Pij-^sA'f) i^ the spectral function of Yiij.^sA'f)- The contribution of the i^(1460) to 
Pus;A o'^ ^h^ LHS is, in the NWA, ^fj^/^^^^-. 5 (q^ — 'm'Kn^Qo)), with the usual Breit-Wigner 
generalization to finite width. The decay constant, /ii-(i46o); is defined as usual by 

(0|A-|X(1460)(g)) = iV2fKiim% ■ (19) 



The analogous NWA contribution of the 7^(1460) to the RHS of Eq. ([TD is 
2/l:(i46o) '^i:(i460) ^ V'?^ ~ ""^/cCueo))- ^ similar relation holds between the 7^0(1430) con- 
tributions to Pulv ^^^ ("^s ~ "^m) Pus;S- It is thus possible, in principle, to determine the 
longitudinal resonance contributions of the last section indirectly, by fixing /x(i460) and 
/i^*(i430) through analyses of the us pseudoscalar and scalar correlators, respectively. In this 
section we will show that such a determination is, indeed, feasible. 

The method involves the analysis of the correlators of the divergences of the us V and 
A currents, U^s-qv = {nXg - niuf n«s;5 and Iius;dA = {rris + rriuf n«s;P using finite energy 
sum rules (FESR) of a type ( "pinch- weighted" ) known to allow an accurate reconstruction 
of the isovector vector spectral function using as input only the OPE, together with PDG 
values for the resonance masses and widths P3|JS^ . We will refer to these correlators as the 
strange scalar (SS) and strange pseudoscalar (SPS) correlators in what follows. 

The basic idea of the analysis is straightforward. Analyticity leads to the general FESR 
relation, 



ac;c;r,c;dspus;dv,dA{s)w{s) = -—. i dsIlus-av,dA{s) w{s) , (20) 

j^bb,bfb 2t:i J\s\=so 



valid for any function w{s) analytic in the region of the contour. The LHS is determined 
by the decay constants of the relevant scalar or pseudoscalar resonances, while the RHS, for 
large enough sq can be evaluated using the OPE. For the case of the analogous isovector 
vector correlator, it has been shown that, although the breakdown of the OPE near the 
timelike real axis for sq ~ m^ is not negligible (so that, for example, FESR's with w{s) = s^ 
are not well-satisfied for sq ~ m"^ H^), even a single zero in w{s) at s = sq is enough 
to produce FESR's that are very well satisfied when the OPE representation for Hud-y is 



employed in the analogue of the RHS of Eq. ( pOD |35[| . A physical understanding of the 



origin of this behavior is provided by the arguments of Ref. [^. As shown in Refs. |l35| , |36 
working simultaneously with FESR's based on the weight families 



WN{y,A) = (1 
WD{y,A) = {l 



2/)(l + ^y) 
y)\l + Ay) 



(21) 



where now y = s/sq, and A is a free parameter, gives enough variability in the weight profile 
to strongly constrain the resonance decay constants in a given channelQ. We will denote the 
class of pinch-weighted FESR's as pFESR's in what follows. 

The OPE representations of the SS and SPS correlators are known up to dimension 
D = 6 P,|5U|, the D = part being determined to four loop order {0{al)) ||2S[. It is 
convenient to work with the second derivative of 11 with respect to Q^, which satisfies a 
homogeneous RG equation, allowing logarithms to be summed by the scale choice fi^ = Q^. 
One has, for the resulting OPE representations @,^, 

3 {rris ± m„)2 ■" 



n" 



us;SPS/SS 



iQ' 



^us;SPS/Ss{Q ) 



^^us;SPS/SS 



iQ' 



D=0 



D=2 



D=A 



i-K^ 



Q' 
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2\3 



1 + -iraiQ ) + 14.1793a(Q^)^ + 77.3683a(g^) 



3 {nis ± muYm 



2 r 



77 



C(3) a(Q 



T2 1 



^ 28 ,^2^ /8557 



23 ,^2^^ / V 1 A 121 , „, 



2\2 



'G 



+ ( 1 + f a(Q" 




n'' 



us;SPS/SS 



iQ' 



{nis ± niuf 



D=6 



Q' 



32 2 



±3 [{rriigqja ■ Gqj + rrijgqia ■ Gqi)] 



(22) 



^This is true even for correlators for which the spectral function contains significant background 
contributions near threshold. As an example, consider the us scalar channel. In Ref. [p^, an ansatz 
for the corresponding spectral function has been constructed, employing the Omnes representa- 
tion for the timelike scalar Kit form factor in combination with certain additional assumptions. 
The resulting spectral function displays a very significant background contribution near threshold 
(associated with the strongly attractive s-wave / = 1/2 Kn interaction) which cannot be well 
represented by the tail of the -fCo(1430) resonance. If one takes as input on the OPE side of the 
wn and wo FESR's the value of rris obtained from a FESR analysis using this spectral ansatz as 
input, and then, with the OPE representation so fixed, makes an incoherent-sum-of-Breit-Wigner- 
resonance ansatz for the spectral function, and uses matching to the OPE sides of the set of wn 
and wd FESR's above to fix the resonance decay constants, one finds that the ^"0(1430) peak 
of the spectral function is reproduced to within ~ 2%, despite the fact that the near-threshold 
region is, of course, not well-reproduced. The reason is obvious: the spectral function is small in 
the near-threshold region and integrals over the spectral function are sensitive dominantly to the 
regions where it is large, i.e., to the regions of the resonance peaks. 



10 



where Iq and Ig are the RG invariant versions of the gluon and strange quark condensate, as 



defined in Ref. [^, pvsA describes the deviation of the four-quark condensates from their 
vacuum saturation values, and the upper (lower) sign corresponds throughout to the SPS 
(SS) case. 

One should bear in mind that, on the theoretical side of SS (SPS) sum rules, the contri- 
bution of direct instantons to the SS (SPS) correlator is not contained in the OPE represen- 
tation. Such effects are known to play a potentially important role in scalar and pseudoscalar 
channels [^0| -|43[|, particularly at lower scales ~ 1 GeV.0 As a result, one must include an 
estimate of direct instanton contributions, in addition to OPE contributions, in the the- 
oretical representation of the SS (SPS) correlator. A convenient, and phenomenologically 
constrained, model for making such an estimate is the instanton liquid model (ILM) |]4^p|- 



ILM contributions to the theoretical side of polynomial-weighted SPS pFESR's can be ob- 



tained from the result 45 



(23) 



where pi ^ (1/0.6 GeV) is the average instanton size (a parameter of the ILM), rj^s is an 
S'C/(3)-breaking factor whose value in the ILM is ~ 0.6 |^2|], and the result is relevant to 
scales ~ 1 GeV^. The corresponding result for the SS channel is obtained by the replacement 
{rris + rriuY -^ -{mg - rriuY. 

It is important to remember that, for a given scale, the ILM contribution to a typical 
scalar or pseudoscalar FESR is much larger than that to the corresponding Borel sum rule 
(BSR)[]. At the scales we will be employing, ILM contributions to the SPS and SS BSR's 
are, in fact, quite small, while those to our pFESR's are still non-negligible. Consistency 
between pFESR and BSR analyses thus represents a non-trivial constraint on the reliability 
of the ILM representation of instanton contributions [^. In Ref. [^, this consistency 



check was implemented as follows. First, the families of pFESR's noted above were used to 
make a simultaneous determination of m^ -|- rrtj and the resonance decay constants relevant 
to the flavor ij pseudoscalar channel, the values obtained for mi + tjij and the resonance 
decay constants being sensitive to whether or not the OPE was supplemented with ILM 



^The Borel transform of the OPE representation of the correlator of the flavor ij pseudoscalar 
density, for example, displays the wrong dependence on the Borel mass, M, in the chiral limit: 
while the exp(— s/M^)-weighted spectral integral must become independent of M in this limit, the 
Borel transformed OPE representation displays a strong dependence on M. 

^The incorrect M dependence of the theoretical side of the Borel transformed pseudoscalar correla- 
tor sum rule is cured once the OPE representation is supplemented with ILM contributions |41,42|. 



The reason is straightforward: ILM contributions to the scalar and pseudoscalar correlators, 

2 

. The modulus of the MacDonald function 



nii;S,p(Q^)) are proportional to Q^ K^i (pi\/Q^] 
i^_i(p/-\/Q^), on the circle Q^ = |Q^|e*^, is typically much larger for non-zero 9 than it is for the 
spacelike point, 9 = 0. The integral around the circle \s\ = sq present on the theoretical side of a 
FESR thus samples regions of the complex Q^-plane where the ILM contributions are enhanced. 
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contributions. The resulting pFESR-generated values for the decay constants were then used 
as input to a BSR analysis of the same correlator, an alternate determination oirrii+mj being 
the output of this analysis. The pFESR and BSR determinations of rrii + Tnj should then be 
consistent. The only non-trivial sensitivity to ILM contributions in the BSR analysis is that 
associated with the input pFESR values for the resonance decay constants. Consistency 
of the two determinations was found only when ILM contributions were included on the 
theoretical sides of the pFESR's [Q. 

In what follows, therefore, we will employ the ILM to estimate direct instanton effects, 
and determine the strange scalar and pseudoscalar resonance decay constants in a combined 
u!n, ujd pFESR analysis. Compatibility of the pFESR and BSR quark mass determinations 
will be imposed as an additional consistency requirement |. The Borel transform of the ILM 
contribution to the SPS correlator required for this consistency check is given by 

3p?K + mO^M^ [iro(p?MV2) + ir,(p?MV2)] . (24) 

That for the SS correlator is obtained by the replacement {rris + rnuY ~^ —{f^s — ^t^uY- 
Expressions for the Borel transform of the OPE representations are well known, and can be 
found in Refs. 0,|29|,|^. 



We use the following input values on the OPE+ILM side of our sum rules: pi 
1/(0.6 GeV) |g,|l, asijnl) = 0.334 ± .022 @,|8l, (a.G^^ = (0.07 ± 0.01) GeV^ [|9 



{mu + md){uu) = -flml (the GMOR relation), 0.7 < {ss) / {uu) = r^ < 1 ||,|9 
{gqaFq) = (o.8 ± 0.2 GeV^) (qq) |5| and < pvsA < 10. The D = 0, 2 and 4 OPE inte- 



grals are evaluated using the contour- improvement prescription p6| , |3T| , since this is known 
to improve convergence and reduce residual scale dependence |^. The running coupling and 
running mass required in this procedure are obtained using the 4-loop-truncated versions of 
the (3 [^ and 7 |^ functions, with the value of ^^(m^) given above as input. 



The analysis of the SPS channel has already been performed in Ref. [^], to which the 
interested reader is referred for a detailed discussion. The results of that analysis are 

m,(2 GeV) = 100 ± 12 MeV (25) 

fKiimo) = 21.4 ±2.8 MeV (26) 

< /i^(i830) < 8.9 MeV , (27) 

where the errors have been obtained by combining the "theory" and "method" errors of 



^Errors associated with uncertainties in the input resonance masses and widths and the input 
values of parameters appearing on the theoretical sides of the sum rules occur for both the pFESR 
and BSR analyses and are strongly correlated. The BSR analysis has additional errors associated 
with the use of the "continuum" approximation for the high-s part of the spectral integral and 
the uncertainty in the choice of the "continuum threshold" parameter. For the SPS case, these 
were estimated to produce an uncertainty of ~ 9% in mg + ruu |Q. We have employed this same 
estimate for the additional BSR uncertainty in our combined pFESR/BSR SS channel analysis. 
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Ref. [0 in quadrature^. The lack of a strong constraint on /ii-(i83o) is a result of the 



smallness of the i^(1830) contribution to the various pFESR spectral integrals. Since only 
the tail of the /^(ISSO) contributes to hadronic r decay, and the endpoint region is strongly 
suppressed by the kinematic weight factor, this uncertainty plays a negligible role for our 
purposes. Two further points should be stressed: first, the value of rris obtained from 
the SPS pFESR/BSR analysis is consistent with that obtained from recent analyses based 
on hadronic r decay data and, second, even if one completely neglects ILM contributions 
(ignoring the resulting inconsistency between pFESR and BSR mass determinations), one 
obtains a value /ii:(i46o) = 22.9 ± 2.7 MeV compatible with that given in Eq. (^) within 
errors]^. For the purposes of determining the strange pseudoscalar longitudinal contributions 
in hadronic r decay, the result of Eq. ( p6|) thus appears very robust. 

A simultaneous pFESR determination of rus + rriu and the excited resonance decay con- 
stants in the SPS channel is possible only because one part of the spectral function (the 
K pole contribution) is well determined experimentally. Unfortunately the experimental 
spectral constraints are considerably weaker in the SS channel. 

The SS spectral function should be dominated by contributions from Kit intermediate 
states up to and including the Kq^USO) region, since the Kq^IASO) displays essentially 
no inelasticity. Unitarity and the Omnes representation of the timelike scalar Kii form 
factor allow one to represent the Kn component of the spectral function in terms of K^^ 
data and Kit phases @J3^. There are, however, ambiguities in this representation. In 
the literature, it has been assumed that a possible polynomial prefactor is absent from the 
Omnes representation and, in addition, that the corresponding asymptotic behavior of the 
Kn phase required by quark counting rules has already been reached at the upper edge of 
the currently accessible experimental range, s ~ 2.9 GeV^. The spectral ansatz which results 
from these assumptions [^ serves as the basis for a number of recent sum rule of analyses 
of the SS channel |jl| and corresponds (reading from Fig. 2 of Ref. [^) to the constraint 



26.2 MeV < /x*(i430) < 31.0 MeV . (28) 

The corresponding value for m^ is, averaging the errors quoted in Refs. ^{d) and p|(b), 

m,(2 GeV) = 115 ± 15 MeV . (29) 

It turns out, however, that quite sizeable deviations from the asymptotic value of the phase 
are allowed in the region of the Kq{1950) without violating the known ChPT constraint on 



^The method errors refer to changes in the output produced by varying the sq and A ranges 
used in the pFESR analysis, or by performing wn or wd family analyses separately, rather than 
a combined analysis. A breakdown of the contributions to the combined error may be found in 



Ref. |46|. 



^'^The biggest impact of neglecting ILM contributions is on fK{i830)j which becomes 14.5ibl.4 MeV. 
The pFESR value of ms{2 GeV) is also altered, the central value becoming 116 MeV, but it is 
difficult to assign meaningful errors to this number since the pFESR and BSR determinations are 
not consistent in this case. 
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the slope of the form factor at s = 0. These can, in turn, produce non-trivial deviations of 



the spectral function from that obtained in Ref. |38], even in the region of the K^{1A30). 
There are thus potentially significant uncertainties not yet reflected in the range of values 
for /i^*(i430) given in Eq. (||). 

Without fully constrained experimental values for the SS spectral function, the pFESR 
and/or BSR analyses allow us to determine fK*{i43o) and fK*{i95o) only after an input value 
for [nis — niuY has been provided on the OPE+ILM side of the sum rules. The reason is 
that, at the scales employed in our analyses, those terms in the OPE proportional to m^ are 
numerically tiny, so the OPE representation is, to a very good approximation, proportional 
to {nis — rriu)'^- Thus, once one finds an optimized spectral ansatz for a particular value 
of TUs — TUu, say rus — rriu = mo, an equally-well-optimized ansatz for any other value, 
nis — f^u = mi, can be obtained simply by rescaling the fitted decay constants by mi /mo- 
The pFESR analysis thus allows only a determination of the ratios fi/{m,s — m,^). 

In our SS pFESR analysis, our spectral ansatz consists of an incoherent sum of i^Q(1430) 
and ^^'0(1950) Breit-Wigner resonance forms, with PDG2000 values of the resonance masses 
and widths. We employ the same pFESR analysis window as used in our earlier study of 
the SPS channel, namely 3.0 GeV^ < sq < 4.0 GeV^ and < A < 4. The different A values 
correspond to weights with significantly different relative weightings between the first and 
second resonance regions, and hence are useful in tightening constraints on the resonance 
decay constants. As noted above, as a self-consistency check on the ILM representation 
of direct instanton effects, we also require consistency between the value of m,s — m,u used 
as input to the pFESR analysis and that obtained as output from the corresponding BSR 
analysis, in which pFESR values of the resonance decay constants are used. The results of 
this determination are 

/....o,Hl7.6±2,0|(!!;f|^). (31) 

where we have combined all sources of error in quadrature. If we consider the range of m,s 
values given in Eq. (pQ]), the corresponding range of fK*(u30) allowed by Eq. (|30| ) is 

20.4 MeV < /i^.(i430) < 32.0 MeV , (32) 

which turns out to be in good agreement with that given by Eq. (|28|). 

With the values above for the decay constants of the SS and SPS resonances it is straight- 
forward to compute the expected r —>■ A'q(1430)z/^ branching fraction, and also the values of 
the resonance contributions to AR^'^^ 
the result of Eq. ( |30D corresponds to 

0.00003 < B{t -^ K*(1430)z/^) < 0.00011 , (33) 

and hence satisfies the constraint given by the PDG2000 upper bound. As expected on 



. Taking 83 MeV < m,(2 GeV) < 130 MeV |T9 
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kinematic grounds, the ii'(1830) and 7^0(1950) contributions to AR^''^^ are negligible]^. 
The corresponding i^(1460) and 7^9(1430) contributions, which follow from Eqs. (^) and 
(0), areQ 



[Ai?(.0'0)j / ' = 0.0052 ± 0.0014 



K{U60) 
L 



and 



A^(o,o) 



i^o*{1430) 



[0.0059 ±0.0011] 



m,(2 GeV) 
100 MeV 



(34) 



(35) 



The sum of the SS and SPS resonance contributions is thus ~ 10% of the K + it pole 
contribution. This level of suppression of resonance relative to pole contributions is in 
fact quite natural, and represents a combination of chiral and kinematic effects. With 



Vk 



mj^/ml. and y^es — (l-4 GeV^j /ml 



, one has w\^ (i/k) = -13 and w}^ {y^ 



.16. 



Thus, although w\^ (yx), which is proportional to mj^, and hence of 0{ms) in the chiral 
counting, is formally suppressed by one power of m^ relative to w}^ (yres), the factor (1 — y)^ 
in w\^ (y) produces a kinematic suppression of w}^ (yres) which largely undoes this effect. 



As a result, one expects resonance contributions to 



3Jres 

Ai?(o,o) 

2 



to be suppressed, relative to 



the K contribution, by the ratio /ii-*{i43o),i^(i460)//K , which has an ml chiral suppression, 
~ 0.1. Naive chiral counting would have produced instead a less strong suppression, of order 
m^, ~ 0.3. 

As noted above, the result of Eq. ( ^l)) is more than an order of magnitude smaller than 
the lower bound implied by the combination of the PDG2000 upper bound on the i^Q(1430) 
branching fraction and the assumption that the OPE representation of the longitudinal (0, 0) 
spectral weight contribution is reliable. To satisfy the lower bound, one would require a value 
of /ft:{i460) a factor of ~ a/20 larger than that given above. Such a large value, however, 
leads to an exceptionally poor "optimized" pFESR OPE/spectral match. The value of irtg 
corresponding to this "optimized" match, moreover, produces a (0 + 1) OPE contribution 
which already exceeds the experimental value for AR^'^\ and hence violates the positivity 
of the longitudinal SS and SPS contributions. 



IV. SUMMARY AND DISCUSSION 

We have shown that determinations of mg based on inclusive (A;, 0) spectral weight anal- 
yses of flavor breaking in hadronic r decay have an unavoidable, unphysical dependence 
on fc, and that the impact of this unphysical behavior on the extracted values for irts is 



^^The ii'o(1950) contribution is a factor of ~ 20 smaller than the -R'o(1430) contribution, the 
iir(1830) contribution a factor of ~ 60 smaller than the iC(1460) contribution. 
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For definiteness, we have computed these contributions using the value \Vus\ = 0.2196 ± 0.0023 
obtained from the analysis of K^^, decay data. The results, of course, scale as l/|V^sp. 
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numerically significant. The problem has been shown to result from the unphysical behavior 
with respect to k of the OPE representation of the longitudinal contributions to the ud and 
us correlators, when truncated at D = 6. If truncation at D = 6 is justified, then the 
problem lies with the D = 2 part of the OPE representation (whose integrated contour- 
improved series, in any case, already displays rather bad non-converging behavior). Part 
of the problem, however, may lie in the neglect of higher dimension contributions, particu- 
larly since contributions un-suppressed by additional factors of as, and having dimensions 
up to D = 8,10 and 12, are in principle present for the (0,0), (1,0) and (2,0) analyses, 
respectively. Since nothing is known, phenomenologically, about the values oi D = 8 and 
higher condensates, the only way to investigate this question is to consider spectral weight 
(or other pFESR) analyses with sq ^ m^, and try to use the sq dependence to separate con- 
tributions of different dimension[^ To work with sq 7^ 171"^, however, one must necessarily 
perform a non-inclusive analysis, since different kinematic factors, both of which are specific 
to So = rn'^, are associated with the 0-1-1 and longitudinal contributions to the experimental 
number distribution. In the region of the spectrum where the separation into + 1 and 
longitudinal contributions is not straightforward (the excited SS and SPS resonance region) 
we have provided determinations of the SS and SPS resonance decay constants accurate to 
~ 10%. This allows one to evaluate the resonance part of the longitudinal contribution with 
an accuracy of ~ 20%. Even for the (0, 0) analysis, such an uncertainty corresponds to only 
a ~ 2% uncertainty on the total longitudinal subtraction; for weights which more strongly 
suppress the excited resonance region, the corresponding uncertainty is even smaller |^. 

We conclude with a comment on the implications of our results for future r decay de- 
terminations of rus. It is our opinion that the bad behavior of the OPE representation of 
the longitudinal contributions precludes the reliable use of an inclusive analysis, and forces 
us to make a theoretical evaluation of the longitudinal resonance contributions to the spec- 
trum. As a result, the previous dis- incentive to studying the sq dependence of any particular 
pFESR (the non-inclusivity of such an analysis) is no longer in play. Since the D = 2,4,6 
terms in the OPE representation of An^°+^) are well-behaved, a study of the sq dependence 
then becomes crucial either to demonstrating explicitly that higher dimension contributions 
can, indeed, be safely neglected or to constraining their magnitude, if they cannot. Because 
of the very strong correlations between spectral integrals corresponding to different Sq, but 
fixed weight, w{s/so), truncated OPE representations which either miss, or pass obliquely 
through the experimental error band for the so-dependent spectral integral results, will both 
signal the presence of such neglected higher D terms. Our determinations of the SS and 
SPS resonance decay constants make it possible for the Sq dependence of the + 1 sum rules 



^^Contributions of dimension D = 2N scale with sq as 1/sq , up to logarithms. 

^^ Although the large amount of new r data that will be generated by the B factory experiments 
will eventually dramatically change the experimental situation, at present, experimental errors on 
the V + A us number distribution are quite large (~ 20 — 30%) beyond the K* region. As a result, 
reduced errors on rris, at least at present, require use of weights which fall off with s in this region 
more strongly than does the (0, 0) transverse weight. The uncertainties on our determinations of 
the decay constants thus play a negligible role in current analyses. 

16 



to be studied in a straightforward manner, and we believe that such a study should be part 
of all future investigations. 
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TABLES 



,(fc,0) 



TABLE L Comparison of longitudinal and total OPE contributions to ARr '"^ as obtained in 
Ref. [10]. The quoted values of ms{m'^) are those obtained in Ref. [10] using the given spectral 
weight, {k,0). This means that the difference between the total longitudinal OPE contribution 



and the experimental value, 



AR 



(fc,o)- 



exp 



, is the value of the + 1 OPE contribution produced by 



the given ms{m'^). The decrease of the longitudinal D = 2 contributions with A; is a reflection of 



the decrease in the extracted m. 



,m: 



with k; as explained in the text, for a fixed ms{m^), the 



longitudinal D = 2 values would be increasing with k. The experimental errors are those of the 
ALEPH collaboration. 



Weight 



ms{ml) [MeV] 



AR 



(fc,o) 



D=2 
L 



AR 



(fc,o) 



OPE 
i 



Ai?('='°) 



exp 



(0,0) 
(1,0) 
(2,0) 



143 
121 
106 



0.201 ±0.085 
0.150 ± 0.073 
0.126 ± 0.072 



0.274 ± 0.087 
0.223 ± 0.076 
0.199 ±0.074 



0.394 ±0.137 
0.383 ± 0.078 
0.373 ± 0.054 



TABLE II. Longitudinal 'k+K pole contributions to ARr 



?(k,0) 



together with the resonance 



contributions implicit in the longitudinal OPE representations of Ref. [10]. The column labelled 



'IND' gives the latter results with each 



ARY 



(fc,o) 



OPE 



evaluated using the central value of the 



corresponding independent fit for ms{rri^). The column labelled 'COMB' gives the same results. 



except that now the 



AR'; 



(kfi) 



obtained in the combined ana^ 



OPE 



are all evaluated using the common value msirn^) = 119 MeV 
ysis using the (0,0), (1,0), and (2,0) spectral weights. 



Weight 



AR. 



(fc,o) 



K+TT 



i. 



^j^m 



resonance 



L. 



(IND) 



AR. 



(fc,o) 



resonance 



L. 



(COMB) 



(0,0) 
(1,0) 
(2,0) 



0.1204 
0.1105 
0.1014 



0.154 ±0.087 
0.112 ±0.076 
0.097 ± 0.074 



0.092 ± 0.062 
0.107 ±0.074 
0.130 ±0.092 
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